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Abstract 



We investigate the eigenvalues statistics of ensembles of normal random matrices when their 
order N tends to infinite. In the model the eigenvalues have uniform density within a region 
determined by a simple analytic polynomial curve. We study the conformal deformations of 
normal random ensembles to Hermitian random ensembles and give sufficient conditions for 
■ the latter to be a Wigner ensemble. 

OO 
t> 

1 Introduction and Statement of Results 

o 

q I Since early fifties Hermitian random matrix theory plays an important role in the statistical de- 
scription of the spectra of complex systems [H |2J. Recently non-Hermitian random matrices have 
been used to treat problems in superconductor physics with columnar defects [31 H] , in quantum 
chaotic systems [3] , and in quantum chromodynamics [HI E] ■ 



Normal random matrix ensembles have been playing a major role in several areas such as in 
the study of fractional quantum Hall effect [5], quantum Hele-Shaw flows |9J, integrable hierarchies 
[TU] , and integrable structure of the Dirichlet boundary problem [TTI [T2] . 

In the present work the ability of normal ensembles to be conformally deformed into Hermitian 
ensembles is exploited to project density of eigenvalues of non-Hermitian matrices into the real 
axis. The problem is addressed by using the so-called invariant ensemble model, characterized by 
the probability of finding a. N x N matrix M of a class within the ensemble given by 

P(M)dM <xex-p{-NTr[V(M)]}dM, (1) 

with the trace Tr[V(M)] and the Riemann volume dM invariant under unitary transformations. The 
corresponding eigenvalue density, in the limit N — > oo, depends on the particular form of V(M). For 
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the Wigner ensemble of Hermitian matrices with V(M) = —M*M (M* is the Hermitian conjugate 

(7 

of M), the entrances of M are independent and identically distributed Gaussian random variables 
with zero mean and variance a 2 /N. The density of eigenvalues follows the Wigner semicircle law 
supported on [— 2a, 2a] [2]: 

d/JL W (x) = sjA- X 2 /a 2 X[-2a,2a] (x) dx, (2) 

where Xa{%) — 1 if x £ A and otherwise. 

A particularly interesting potential has been put forward by Wiegmann, Zabrodin and coworkers 
[TUl [TlT [T2| [T3] who established a connection between normal random matrices and conformal 
mappings. They considered 

V(M) = —(M*M - p{M) - p{M)*), (3) 

where 

p(z) = Y J t^ (4) 

3>l 

with t > and tj G C. As N — > oo, they showed, at the level of formal manipulations, that : (A) 
the density of eigenvalues is uniform within a simply connected domain D C C whose boundary is 
given by a simple analytic curve 7; (S) the domain D is characterized by the fact that its exterior 
harmonic moments 



tj = (t zz J dz, j > 1 , (5) 

where 7rt stands for the area of D, are the coefficients of (j3J); and (C) the Riemann mapping from 
the exterior of the unit disk onto the exterior of the domain D obeys, as a function of the tj, the 
equations of the integrable dispersionless Toda hierarchy. 

Potentials of the form (J2D give rise two sorts of mathematical problems. Except in the case of 
polynomial p(z) of degree 2, where the domain D is bounded by an ellipse, V(z) is not bounded 
from bellow and integrals with respect to (CQ) diverge. The other problem concerns with the fact that 
D may not be uniquely determined by the moments (jSJ). From the point of view of equilibrium 
measures (see Section [2]), a relevant fraction of eigenvalues of a M may escape to infinity or to 
another Riemann surface. 

Recently, the results (A) and (B) have been set in a rigorous frame by Elbau and Felder [14] . 
To avoid the above mentioned problems, they consider the following restrictions: 

Elbau-Felder Potential. If 

p(z) = t lZ + t 2 z 2 + ■■■ + t n+1 z n+1 (6) 

is an analytic polynomial of degree n+1 with t > and t = (ti, . . . , t n +i) £ C n+1 such that t\ = 0, 
1^2 1 < 1/2, Elbau-Felder potential is a real-valued function on C given by 
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V{z) = y{\z\ 2 -p{z)-p{z)) . 
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It can be shown by direct computations that V(z), under the above conditions, is positive in a 
neighborhood of z = and has a non-degenerate absolute minimum at z — 0. From now on, V{z) 
shall stand for the Elbau-Felder potential. The problem of divergence of the integrals is solved 
by Elbau-Felder in a naive way - imposing that the eigenvalues of matrices within the normal 
ensemble remains bounded: 

Elbau-Felder Ensemble. Let £ C C be the closure of a bounded open set that contains the 
origin and consider the following class of matrices 

M N {T) = {Ae Mat c (iV) : [A, A*} = 0, a (A) C £} (7) 

where cr(A) denotes the spectrum of A. An ensemble is said to be of Elbau-Felder type of degree 
n + 1 if it fulfills conditions stated between and ©. A closed polynomial curve 7 of degree n 
can be parametrized by 

n 

w 1 — > h(w) = rw + ajW~i , \w\ — 1 (8) 

j=0 

for some r > and the aj G C. Elbau and Felder have shown that, as long as < 1/2 and 
to is small enough, the problem of determining the exterior moments tj out of the curve has a 
unique solution for simple closed analytic polynomial curves. They give a set of equations that 
defines an invertible map F : (r 2 , a , . . . , a n /r n ) — ► (t , . . . , t n+ i) from R x C n+1 into itself about 
(0, 0, 2i 2 , . . . ,{n + l)i n+ i) {tj stand for the complex conjugate of the tj). By the Euler-Lagrange 
variational equations, the eigenvalues density is uniform in D. We refer to Theorem [2] for a precise 
statement. 

In the present work, we study conformal deformations of the Elbau-Felder ensembles into Her- 
mitian ensembles. This is achieved by considering a family of polynomial curves of degree n: 
w 1— > h{w; s), with the a,j(s) depending on a parameter s G (0, 1]. This family is chosen in such 
a way that h(w; 1) = h(w) parametrizes the initial curve 7 whose interior domain D supports the 
eigenvalues. After the construction of h{w; s) the support D (resp. the harmonic moments tj) also 
depends on s under s 1— > D(s) (resp. t 1— > tj(s)). To state our result, we denote by r = (ti, . . . , T n+ %) 
a vector on the affine space Z C C n+1 with n = and \t%\ = 1- 

Theorem 1 Consider the Elbau-Felder ensemble with to > and s 1— > t(s) G C n+1 snc/i that t\ = 0, 

*2 = — 2 — exp ^ s ' 

and 

t j = s Aj Tj for 3 < j < n + 1 , 

wift <y9 G [0,27r) ; Tj G C, Aj > 1 and s G (0, 1]. There exist r = r (r) > such that for every 
< r < r : 
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1) There is a unique simple analytic closed polynomial curve 7 = j(s,r,T) of degree n, with 
external harmonic moments t(s) and area of interior domain irt with t depending on t, r 
and s. 

2) The curve is parametrized by h(w; s,r, r) = rw + YTj=o r ^ a j w ~^ w ^ \ w \ = 1 an ^ a j = 
ctj(s, t) is uniquely determined by t and t(s). 

3) The eigenvalue density is uniform within D, the interior domain of j, for every s £ (0, 1]. 
Moreover, if Aj > 1 then the Elbau-Felder ensemble can be conformally deformed, as s goes 
to ; into a Wigner ensemble with support on [— 2r, 2r\. 

Remark 1 As long as < r < r , h(w; s,r, r) is a Riemann mapping from the exterior of the 
unit disk onto the exterior domain D of 7 and the area irt of the domain D remains positive for 
alls £ (0,1]. 

Remark 2 The assumption |i 2 | < 1/2 in Theorem [H breaks down when the exterior domain 
D_(s) = C\D(s) is deformed into the slit domain C\[— 2r, 2r]. In Section^we generalize Elbau- 
Felder's results using Crandall-Rabinowitz bifurcation theory from simple eigenvalues (see e.g. fTh^ ) 
to construct a parametrization that allow us to lett2 — > 1/2 maintaining the parameter r away from 
0. Elbau-Felder 7 s parametrization, coming from the implicit function theorem applied to the map 
F , mentioned right below Eq. |2P, defines a curve in R x C n+1 that bifurcates at ti = 1/2. 

Remark 3 If Aj > 1 then the aj = atj(s,T) behave, for s —> 0, as 

(1) a (s) = o(s), 

(a) = 1 — s/2, 

(Hi) aij(s) = o(s) for 1 < j < n. 

This relations will be used to prove item 3) of Theorem^ 

Remark 4 Mashkov et al (Sec. 6 of JU| /j considered a family of simple closed analytic curves 7(5) 
given by an equation P(x,y/s) = 0, w = x + iy, converging to the segment of the real line [a, f3] as 
s —>■ 0. At the level of formal manipulations, they have shown that the eigenvalue density projected 
into [a, 0\ yields 

p{x) = lim Ay ^' J ^ = y/(x-a)(P-z)M(x) 

where Ay(x; s) is the width of D(s) at the x coordinate and M(x) is a smooth function, regular at 
the edges. 
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This paper is organized as follows. Section [2] presents some preliminary results and introduces 
the two ingredients, the balayage problem and the Schwarz function, required for the proof of 
Theorem [TJ In Section [3] we prove two auxiliary results, Propositions [1] and [2j Section H] uses 
Crandall-Rabinowitz bifurcation theory from simple eigenvalues to establish a smooth inverse map 
F^ 1 in 1 x C n+1 about t\ = and = 1/2. Theorem H] in Section [5] gives an explicit expression 
of the Balayage measure for the potential V. Section [6] concludes the proof of Theorem [1] based in 
Lemma [1] and Section [7] gives some examples. We present in Section [8] our conclusions and Lemma 
[2] is proved in Appendix [A] 



2 Basic Setting 



2.1 Eigenvalue Distribution for Normal Ensembles 

For normal unitarily invariant ensembles, we can write Eq. (pQ) in terms of the spectral coordinates. 
The joint probability of the eigenvalues {zi\f =l C E of M reads 

TV 

P N (zi,--- ,z N ) ocexp{-(2 Yl lo s\zi-Zj\' l + Nj2v(zi))}- (9) 

l<i<j<N i=l 

Introducing the empirical measure of the eigenvalues 

N 

dfi N (z) = N~ l J2d(z - Zi)d 2 z, (10) 

i=l 

([ID can be written as 

^n{zi, ■ ■ ■ , zn) — ^ N e yr 
where Zn is the normalization and 



I» = J (V(z) + U»(z))d»(z) (11) 
is the total energy. The logarithmic potential associated with /i given by 

U»(z)= I \og\z-w\' 1 dn{w) . (12) 



The integrals with respect to (Q have, in the limit iV — >• oo, dominant contribution governed by a 
variational problem: 

E V = inf I V (»), (13) 

where the infimum is taken over the set .M(E) of Borel probability measures in E C C. If a 
probability measure [i v satisfying 

E v = I^ v ) 
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exists, it is called the equilibrium measure associated with V. 
to converge weakly to a unique equilibrium measure as iV - 
and [9] for normal ensembles). 



The empirical measure (fTOj) is known 
■ oo (see [2] for Hermitian ensembles 



Theorem 2 (Elbau-Felder) Consider the Elbau-Felder ensemble of degree n. There is 5 > 
such that for all < t < 5 a unique equilibrium measure d\x exists and is uniform within a domain 
D C S that contains the origin: 

dfi = —Xd (z) d 2 z; (14) 
7rt 

D is uniquely determined by the exterior harmonic moments (ti = 0) 

nto = d 2 z the area ofD 
Jd 

tk = — - [ z~ h d 2 z , if k = 2, • • • , n + 1, 
J C \ D 

t k = Oifk>n + l (15) 

and its boundary 7 is a simple closed analytic polynomial curve of degree n; if h(w) = rw + + 
ai/w + ■ ■ ■ + a n /w n , \w\ = 1, parametrizes 7, then 

n 

t =r 2 -J2jWj\ 2 . (16) 
i=i 

There exist homogeneous universal polynomials Pjk G Z[r, ciq, • • • , a^-j] of degree k — j + 1, 1 < j < 
k < n + 1 swc/i that 

n 

jtj = dj-xr'^ 1 + ^2 a-kr~ k Pjk{r, a , ■ • ■ , a k -j) (17) 

k=j 

is an invertible transformation from RxC" into itself in a neighborhood of(r 2 , ao, a±/r . . . , a n /r n ) = 
(0, 0, 2i 2 , . . . , (n + l)t n+ i). For sufficiently small r, the function h(w) is a Riemann mapping from 
the exterior of the unit disk onto the exterior of D. 



The proof of existence of the equilibrium measure requires to verify Euler-Lagrange type equa- 
tions. Let ^ 

d 2 ( (18) 

be a function defined in S given by V plus the logarithmic potential (1121) associated with the 
uniform measure \x in D. Lemma 6.3 of [14J shows that E(z) = holds for almost every z G D. 
According to Corollary 3.5 of [H], [i is the unique equilibrium measure if 

E{z) > for every z G T\D . (19) 

We extend Elbau-Felder's proof to the near-slit-domains in Section HI 



E(z) = V(z) 
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2.2 Balayage Problem 



We tackle the problem of analyzing conformal deformations of normal ensemble by means of bal- 
ayage techniques [IB] . This allows us to solve the problem focusing only on the behavior of the 
boundary 7 of D. Therefore, in our approach the balayage technique plays a major role. Let G C C 
be an open set and dG its boundary. 

Let v be a probability measure on G (such that v (C \G) = 0) and let the logarithmic potential 
U u (see Eq. ([12]) ) be finite and continuous on G. The balayage problem (or "sweeping out" 
problem) consists in finding a probability measure v with support on dG such that 

U" = U v almost everywhere on dG. (20) 

We call v the balayage measure associated with v. Throughout this work, we consider the following 
space of functions: 

Definition 1 Let G C C be a bounded open set. We denote by 7i(G) the space of all holomorphic 
functions on G and continuous on its closure G. 

If G C C is a bounded open set and v is a probability measure with compact support in G, 
then v is the unique measure supported in dG satisfying (1201 and such that U v (z) is bounded in 
dG. In addition, v possesses the following property (see Theorem II — 4.1 of [PT]): 

[ fdu= [ fdu (21) 

JG JdG 

holds for every / G TC(G). 

We may choose G the interior D = D\y of the compact support D of the equilibrium measure 
\i associated with V. In this case, we write fj, — fi\ff + yu| 7 and sweep out only the part lying on 
G: fl = fi\^ + /i| 7 . Since the equilibrium measure has no mass concentrated in 7, we have /i| 7 = 0. 
The balayage measure associated with the equilibrium measure /1 is denoted simply by /I and has 
support in 7. 

2.3 Parametric Curves and Schwarz Function 

The following definitions will be important for the characterization of the curves appearing in our 
main result. We shall start with the basic 

Definition 2 A curve T in C is said to be simple if there exist a parametrization t 1— > h(t) for 
t G [a, b] such that h(t) is injective, i. e., if for all x, y G [a, b] h(x) 7^ h(y) when x 7^ y. If 
h(a) = h(b), in this case T is said to be a simple closed curve. A curve T in C is said to be an 
analytic curve if there exist a parametrization t 1— > h(t) for t G [a, b] such that h is analytic and 
h'(t) 7^ forte [a,b]. 
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Next we introduce the polynomial curves on the complex plane. 

Definition 3 A curve V in C is said to be a polynomial curve of degree n if it is parametrically 
represented as 

h{w) = rw + a + aiw~ l + . . . + a n w~ n . (22) 
with r > 0, a n ^ and \w\ = 1. 

We shall define the Schwarz function 

Definition 4 Let T in C be an analytic arc and let Q be a strip-like neighborhood of V. The 
Schwarz function S is the unique analytic function on Q such that 

S(z) = z , z £ r. 

For a treatise on the Schwarz function with applications see [THJ \W\ . 

Remark 5 Hereafter, 7 denotes a simple closed analytic polynomial curve. Moreover, S stands 
for the Schwarz function 0/7. 

Schwarz function S will play a major role in the conformal deformation of the Elbau-Felder 
ensemble. We shall show the balayage measure is proportional to the Schwarz function S. 

3 Riemann Map 

We shall prove some auxiliary results, which concern the behavior of the family h(w; s) as the 7(5) 
is deformed to the real line. 

n 

Proposition 1 Let h(w; s) = rw -\- s ^^aj(s)w~\ \w\ = 1, be for each s £ (0, 1] the parametrization 

j=0 

of a closed polynomial curve 7(3) of degree n with h(-; 1) = 7 a simple curve. Let the condition 

n 

£(s):=r-J2j M*)\ >0 (23) 

be satisfied for every s £ (0, 1]. Then, for each s £ (0, 1], the 7(5) remains a simple polynomial 
curves and h(w; s), as a map from the exterior of the unit disk into the exterior 0/7(5), is biholo- 
morphic (a Riemann map). Furthermore, for every s £ [S, 1], < 6 < 1, t (s) is bounded away 
from zero. 



S 



Proof. Let us begin with the estimation of to- It follows from (|23|) that r > |%(s)| holds for every 
j. Multiplying £(s) by r it yields that 

n n 

< r 2 - 3 \ a A s )\ r <r 2 -^] \a 3 {s)\ 2 = t {s) 
j=i 3=1 

is bounded away from zero for s G [5, 1], < 5 < 1. Eq. ( 1231) also implies that /i(u>; s) is an analytic 
curve, that is, the derivative of h(w; s) with respect to w, denoted by h'{w\ s) is bounded away 
from zero: 



\h'(w;s)\ 



(lj{S)W 



-O'+i) 



> r- J^j 



-a+i) 



where in the last passage we have used the triangular inequality Since we are analyzing the exterior 
of the unity circle \w\ > 1, we have 



\h'(w;s)\>r-^2j\a j (s)\>0yse(0,l], 
i=i 

which also holds in a small neighborhood of \w\ = 1. 

Now, for every w, z G C with \w\ = \z\ = 1, by the triangular inequality, 



(24) 



\h(w; s) — h(z; s)\ > r \w — z\ — \J \a>j(s) 

3=1 



1 1 

w 3 z 3 



r \w — z\ — |a_j(s) | \w 3 — z 3 \ 

3=1 



> £(s) \w — z\ > 



(25) 



if w 7^ z. The last passage follows from \w 3 ' — z 3 \ < j \w — z\ which can be shown using the 
telescopic identity 

w 3 — z 3 = w 3 ^ 1 (w — z) + w 3 ~ 2 {w — z)z + • • - + (w — z)z 1 ~ 1 

together with the triangular inequality Equations ff25l) and fT2"4"|) imply that the map h(-; s) : S 1 — ► 
C is an embedding, 7(5) is a simple curve and h(w; s) is a Riemann map from the exterior of the 
unit circle onto the exterior of 7(s) for every s G (0, 1]. The polynomial curve 'j(s) with < s < 1 
preserves all properties assumed for 7(1) = 7, concluding the proof Proposition [H □ 



Remark 6 Proposition^ has an intuitive appeal. To a polynomial curve 7(5) fail to be simple it 
has to develop a cusp. However, when h(w c , s) form a cusp, we have h'(w c , s) = 0, this situation is 
prevented as long as £(s) > 0. Proposition^ gives a sufficient condition for h'(w, s) 7^ and show 
that is also sufficient for 7(3) to remain a simple curve. 
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The next result concerns the conditions (i — Hi) of Remark [3] and the deformation 7(3) to the 
real line as s — > 0. 

Proposition 2 Consider a polynomial curve 7(3) of degree n parametrized by h(w; s) satisfying 
conditions (i — Hi). Then, lim^o h(w; s) maps the exterior of the unit disk onto C\[— 2r, 2r]. 

Proof. As h(w, s) is a Riemann map for all s G (0, 1] by Proposition [I], then it suffices to show 
that lim s _^o7( s ) = [— 2r, 2r]. Indeed, since \w\ = 1 we may choose w = e tB with 6 G [0, 2tt]. The 
Riemann map reads 

■n 

h(e ie ; s) = re ie + a (s) + ai(s)e" ie + ^ a 3 (s)e~ i] \ 

i=2 

which, under the hypotheses, yields 

h(w;s) = r{e i9 + e~ w ) + o{s), 
implying that lim s ^ h(e l9 , s) = 2rcos9 G [— 2r, 2r] for all 9 G [0, 2n]. □ 



4 Exterior Harmonic Moments of Near-to-Slit Domains 

Let t = (ti, t2, ■ ■ ■ , t n+ \) be the exterior harmonic moment of the domain D - containing the origin 
and bounded by 7 - and let irt be the of D. 

When a given collection t of n + 1 complex numbers, together with t > 0, determines a simple 
polynomial curve 7 of degree n? We refer to Theorem 5.3 of [H] for the solution to this moment 
problem. If t\ = and |t2| < 1/2, then every complex numbers £2, ■ ■ ■ ,t n +i determine a curve 
7 with these properties provided to is sufficiently small. We shall give a proof of this result in a 
language more appropriated for the generalization needed. 

The map (p,a) G R+ x C n+1 1 — ► (to,*) G M+ x C n+1 defined by (DED and by the contour 
integral ([5]), 



jtj = -— ® /i(u> 1 )h'(w)h j (w)dw 
2tu ./ h= i 



taking residues at infinity: 



E ^ Res ( ^ \, , v/,„W ; 00 ) (26) 



v 1 + TJi=o^/w i 



tj = _ a i a o - ^1 + + O (p 2 ) , 1 < j < n + 1 (27) 
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with p = r 2 , a,j = r ia,j, < j < n, and aj = if j > n, has a smooth inverse in a neighborhood 
of (0,0, t 2 , ■ ■ - ,t n+ i) provided \t 2 \ ^ 1/2. 

In this section, we show how the inverse function theorem is applied in this situation and extend 
it for the case \t 2 \ = 1/2. We also verify whether the inverse determines a simple polynomial curve 
7 = 3D and whether a measure p, uniform in a near-to-slit domain D, is the equilibrium measure 
of the Elbau-Felder ensemble. 



4.1 Bifurcating Curves 



We observe that p* = = and a* = (j + l)t,+i, j = 1, . . . ,n solve the equations (12?]) for 
(p, ao, . . . , a n ) and because it takes complex conjugation of the tj + x, we look ( 1271) as a map from 
R + x C n+1 x C n+1 into itself 

(28) 

and we write 







F : (j),a,a) i — ► (to, *» *) 
F(p, a* + pep, a* + pep) - (t , t, t) 



(29) 



= KhViv) +p(p, i Pi c p) - (to, o,o) . 

Here, F has been expanded in Taylor series about (p*, ex.*, 6c*) with remainder pp(p, <p, cp) = O (p 2 
r = to/p = O (1) an d / is the linear map 





( 1 ^ 


( 1 


-4|t 2 | 2 


T T \ 


( 1 ^ 


l(l,<p,<p)=L\p*,a*,a*] 






—v 


-K J- 1 






\<P J 


V 


—v 


J- 1 -K j 





(30) 



where is the null column vector in C n , Vj = 2 (1 + (j + 2)f 2 ^+2 if 1 < j < n with 
v n = v n+ i = 0,0 denotes the (n+ 1) x (n+ 1) null matrix, J = diag {j}™=i and K is a (n+ 1) x (n+ 1) 
matrix with = (i + l)tj+i for 1 < i < n and otherwise. 

Since L is invertible for |t 2 | ^ 1/2, (1281) has a smooth inverse defined in a neighborhood of 
(t Q ,t) = (0, 0, t 2 , . . . , t n+ i). The implicit function theorem applied to (T2"9"l) (with t e C n+1 fixed) 
uniquely defines two smooth curves parametrized by p: 

<p(p) = T« + E?; + 0(p) 

r (p) = 1 -4|t 2 | 2 + 0(p) 

on C n+1 and K+, respectively, where B = (1 - 4 |t 2 | 2 ) _1 Jif and T = J + 2t 2J B. We note that the 
leading order in p of equation (|2"9~I) can be written as 

|2 



-if J" 1 

J- 1 -K 



1 - 4 \t 7 



To 



V 
V 



(31) 
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whose solution is given in Appendix lAl 

The function t — P T o{p) is monotone increasing in p G [0, a] for |t 2 | < 1/2 and sufficiently small 
a, and its inverse, p(to, t), is a well defined function of to and t. The inverse of (1271) in R + x C n+1 
thus reads 

(t ,t) ^ (p(t ,t),(x* + p(t ,t) <pop(t ,t)) . (32) 

The above application of the implicit function theorem breaks down if the second harmonic 
moment t 2 tends to 1/2 and this is the case when the external domain of y(s) tends, as s — > 0, 
to the slit domain C/ [— 2r, 2r] (see Proposition [2]). We need, therefore, to extend Theorem 5.3 of 
[Hj to include this case. For this, we shall rescale all components of t, excepted t 2 whose modulus 
square will be denoted by A = |t 2 | , as 

t J = (l-4A)r, ,j^2 (33) 

and apply the constructive bifurcation theory from a simple eigenvalue developed by Crandall 
and Rabinowitz (see e.g. [15]). The theory applied to equation (1291) uses the implicit function 
theorem with the role of p and A exchanged. Instead of the two parametric curves (fj = <fj(p) and 
r o = r o(p), we consider <pj = <fj(p, A) and f = f (p, A) as a function of A for (p, r) G C n+1 fixed, 
where r = (ti, t 2 , . . . , T n +i) is a vector in the affine space of C n+1 with |-7~2 1 = 1, denoted by Z. 

We observe that (p*, a*) with |t 2 | = 1/2 is a bifurcation point for ff271) since every neighborhood 
of this point contains a solution which differs from (132]) . Note that the tangent map of (|27|) at 
(p*,a*) with \t 2 \ = 1/2 is singular i. e., L is not invertible at A = \t 2 \ 2 = 1/4. Using bifurcation 
theory, we shall construct a pair of smooth curves for fixed (p, t) G C n+1 : <p 3 - = <Pj(A) and 
f = r (A), A G (1/4 - 6,1/4] for some b > 0, such that <Pj(l/4) = ^(p = 0, |t 2 | = 1/2) and 
fo(l/4)=T (p = 0,|t 2 | = l/2). 

Proposition 3 Given (p, t) G M + x Z ~ C™ +1 ; t/iere exzst two uniquely defined smooth curves, 
A i — > <^(p, A, t) on C n+1 and A i— > f (p, A, r) on R, defined by §EBj) , such that the inverse of (FJ\ ) in 
a neighborhood of (0, 0, r) zn R + x R + x Z is written as 

(to, P, t) i — > (p, a* o A(t , P, t) + p(a + C ° A(t , P, t))) (34) 

where A (to, P, T ) a function of t , p and r which is the unique solution for A oft — pf (p, A, r) 
m t/ie domain < A < 1/4, p G [0,f 2 ] wrat/i f = r(r) sufficiently small. Moreover, there exist 
f = f(r) > such that h{w), defined by (TJJ] with r = yfp and aj = raj given by (3$) , parametrizes 
a simple polynomial curve of order n which can be deformed to the segment [— 2r, 2r], for every 
fixed r < f. 

Proof. It suffices, for the first statement, to verify the hypothesis of Theorem 1 in Sec. 3 of [15] . 
For this, it is convenient to write ( 1291) as 

/ (1, if, ¥>) + (!- 4A)/x(l, <p, <p) + p(A, p, <p, <p) - (r , 0, 0) (35) 
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where Iq and l x are linear maps in M + x C™ +1 x C n+1 , with l 





/ 1 ^ 


( 










\ <p ) 


\ 



!| A=1/4 and 4/x = - d(l + p)/dX\ x=1/4 : 

( T T \ 

-K J- 1 



( 1 ^ 
\ <P ) 



( 1 ^ 






i 


T 


T 


\ 




( 1 ^ 








-«! 


-K x 











V ^ y 




V 


-Vi 





-K x 


> 







(36) 



h{l,(f, (f) = L x if = -v x -K x O cp (37) 

\ <P ) \-v x O -K x J \ Cp ) 
with (v x ) j = (1 + 1/j) (j + 2)r j+2 /r 2 + 0{p), (Kq)- = r 2 8 ix 8 jX + 0(p) and, given 

w x = (-t 2 /2, 3r 3 , ...,(n+ l)r„,+i, 0) , 

K\ is, up to the leading order in p, a in + 1) x in + 1) matrix with u>i in its first column and 
everywhere else. Moreover, p(A, p, cp, cp) is a smooth map from M + x M + x C n+1 x C rt+1 to 
E. x C n+1 x C n+1 which satisfies 



p(l/4, p, cp) = p(A, 0, v?, cp) = ^(1/4, p,(p,(p) = 0. 

We observe that and consequently fl33|) . extends by continuity to p = 0. In this way, we 
define 

F(p, a* + p(a + C), a* + p(a + 0) - (*o, *) 



G(p,A,C,C) 



for (A, p, C, C) e [I/ 4 - 1/4 + 6] x [0, a] x C n+1 x C n+1 where, for a 



-n+1 v (Pn+l „w w .v. - 0, . . . , 0), (1, a , a ) 

is an eigenvector of Lq associated with the (simple) null eigenvalue. It follows by f l35|) - fl37|) . that 



G(p, A, C, C) = k(0, C, + (1 - 4A)Z 1 (1 ) a + C, <*o + C) + p(A, p, c*o + C, «o + C) - (r , 0, 0) . 

and the implicit function theory can be applied to G = provided the derivative of G(p, A, £) 
about (A, £, £) = (1/4, 0, 0), defined by the linear map 



/ 



(A, C, C) i— io(0, C, - 4^(1, a , a„) 



O 2 



\ 4(«! + wx/y/rfi J- 1 -K Q J 



( X\ 

C 



is nonsingular (see Theorem 1 in Sec. 3 of [H]). Since it is always invertible, G = (with p > 
and r G Z fixed) uniquely defines two smooth parametric curves: 



C(A) = f 2 ^ («i + wi/v 7 ^) +^0(^1 + ^1^) +0(l-4A,p) 
f (A) = l-4A + 0((l-4A) 2 ,p) 



(3? 
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on C n+1 and on M + . As in the previous case, writing 77 = 1 — 4A, the leading order of equation 
G = reads 

V = TO 

-*"7*' J " = ,( (39) 

from where (}38l) can be readily obtained. 

The function to = P r o{p,^) is monotone decreasing in A, for < A < 1/4, p e [0,f 2 ] with f 
sufficiently small, and its inverse, X(t ,p, r), is a well defined function of t ,p and r. The inverse 
of (TJTj) in R + x M + x Z is thus given bv By (j4"Ij) . f = f(r) is defined by 

n 

2E^ + i) 2 ki+il 2 ^- 1, = i. 

To conclude the proof, it remains to verify that 

Ww)=ra + a + r 1 h r n — , liul = 1 (40) 

W7 W n 

with 

a i = M*o, P, t) + p(a ,j + Q A(to, p, t)) (41) 
parametrizes a simple polynomial curve 7 = 7(^0) • For this, in order 7 to approach the segment 
I = [— 2r, 2r], as the area of its interior nto tends to 0, we set r 2 = 1 and let r 3 and p = r 2 so that 

«o = P Ki + Ci) = 123?e(r 3 )r 2 + 0(r 4 , t„) (42) 

remains inside J, closed to the origin. Now, it is enough to verify the hypothesis ff23l) of Proposition 
[TJ It follows immediately from 

«i = ( 1 -^)+°(^o) 



% = ^(j + l)f J+1 + O(t 2 ) , 2<j<n (43) 



for every r, that 



3=1 V j=2 



holds for r < f where f is defined by equating the expression between parenthesis to 0. □ 
Remark 7 r in Theorem [1\ is thus given by min (f , f) . 

Remark 8 The results of Proposition can be adapted for $33\) scaled with different power of 
(1 - 4A): tj = (1 - 4A) A ^r i; Aj > 1 for each j ^ 2. 

Note that Proposition [3] holds for simple closed analytic curves provided I, Iq and l\ are Frechet 
derivatives of F with respect to an appropriate Banach space. 
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4.2 Equilibrium Measure 



Equation (U9l) is equivalent to 

< £(w) := E(h(w)) 



\h(w)\ z - |/i(l)r + 2Sfte / h^mOdt 



(44) 



for w E h 1 (£) such that |iu| > 1 (see eq. (16) of [Ej), and it suffices to verify only for |w| > 1/R, 
where R = max{|w| : h'{w) = , w G C} is the critical radius of 7, and r sufficiently small. 

Let h^(w) = rw + ao + rai/w be the parametrization of an ellipse that approximates 7 and let 
(^o°\ 4°' ') denote the corresponding external harmonic moment. We denote by S^(w) the function 
defined in ( 144"]) for the ellipse and note that, by Subsection 6.2 of [H], t$ £(°)(w) remains bounded 
for w G h (E) such that \w\ > 1. For (1 — to/r 2 )^ 1 < \w\ < r~ a , < a < 1/3, we have 

1 r 2 a 2 



- (h(w' r ) - /^(tiT 1 )) h'(w) 
to 

-h^iw- 1 ) (ti(w)-hW(w)) 
to 



-r + O(t ) = 0(r L - 2a ) 



-2r 2 a 2 



to — 



to w 

1 

— roLiiv- 
to 

^a 2 2 + O(t ) = O(l) 



+ O(t ) = 0(r) 



t £ {0 \w) = 0{r 
uniformly in t as t — > 0. Consequently, 

\S(w)\ > -\t £ {0 \w)\ - \£(w)-£ {0 \w 



2-a\ 



> I|t £(°)( w )|-I \h(w)\ 2 -\hV\w)\ 2 -\h(l)\ 2 + \hW(l)\ 
to tQ 



2 

T Q 
2 



Ke 



WrV^'t 1 )) h'(Qdc 



HC 1 ) (h'(C) - h {0) '(0) dC 



f / (0) 

l o - l 4. c(0)( 



-to£ w (w) 



is strictly positive for i sufficiently small. For w G h 1 (S) such that \w\ > r "we may proceed 
exactly as in |14j . 



5 Explicit Balayage Measures 

To explicitly obtain the balayage measure in terms of the potentials we need an auxiliary result. 
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Theorem 3 If V : E — > K zs a potential defined in a compact set £ C C loii/i continuous sec- 
ond partial derivatives in its interior, then the variational problem HT3\) is attained at a unique 
equilibrium measure fi v supported in a compact set DcS given by 

d^ v (z) = ^-AVd 2 z, (45) 

at almost every (w.r.t. the Lebesgue measure d 2 z) interior point of D , where A = d 2 j dx 2 + d 2 j dy 2 
is the Laplace operator. 

Proof. Using the smoothness of V the proof follows from Theorems 7 — 1.3 pp. 27 and II — 1.3 pp. 
85 of [17J. □ 

We derive an explicit equation for the balayage measure in terms of the potential V. The result 
will be of fundamental importance to establish the conformal deformation of the Elbau-Felder 
ensemble. 

A direct application of this theorem to the potential V shows that the density of eigenvalues 
is indeed uniform within D. The Elbau-Felder potential reads V(z) = [zz — p(z) — p(z))/to, 
remembering that A = 4d z dg and p(z) is analytic (dgp(z) = d z p(z) =0), it follows that dfi(z) = 
l/(7rto)d 2 z within D. Moreover, we have the following result 

Theorem 4 Let V be an Elbau-Felder potential. The balayage measure fx associated with the 
equilibrium measure ii v with support on 7 = 3D is 

d(2(z) = -L-S(z)dz, (46) 
Zmto 

where dz is the measure of arc length on 7. 

Proof. The main ingredient is the Green theorem. Our hypotheses on V, D and dD guarantee 
that the Green theorem is applicable. Using that A = 4d zz - and / G 7i(D), by (|45p and applying 
the Green theorem we have 

f(z)dfi v (z) = - [ c\{f{z)d z V{z))d 2 z 



L> 71 JD 

1 



f{z)d z V(z)dz 

dD 

and by the balayage measure property Eq. ( |2TT) . we have 

f(z)dfl v (z) = -!- / f(z)d z V(z)dz 

dD lm JdD 

holds for every continuous function on 3D, from where it follows the continuity of the balayage 
measure with respect to the Lebesgue measure. Next, we have that 
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the term of the sum does not contribute to a contour integral with respect to ju for test functions 
/ G 'H(D) (by continuity and by the Cauchy theorem). Then, by using the definition z = S(z), 
z G 7, of the Schwarz function we obtain Eq. (j46p . concluding the proof. □ 

Eq. (j4"6l) relates the equilibrium measure with the Schwarz function of boundary of the support. 
We shall conclude, by applying the Cauchy theorem, that only the branch cut in the interior domain 
D of the Schwarz function will contribute to the balayage measure. Thus, questions concerning 
the equilibrium measure \x turns to the behavior of the Schwarz function. It turns out that, except 
when 7 is a line or a circle arc, the Schwarz function S always has a branch cut [18J. Our next 
result draws some conclusions on the behavior of the branch cuts. 

Proposition 4 If 7 is a simple closed analytic curve, then the Schwarz function S associated with 
7 must have branching points in its interior. 

Proof. Take / G T~i(D) such that f D f(z)d[i 7^ 0. By the property (|2Tj) of balayage measure together 
with the Theorem HJ we have 



Suppose now that S has no branch point in the interior of 7, then by using the Cauchy theorem 



have an even number of branch points inside D. Another important result is that the branch point 
never touches the curve 7. This can be proved by noting that 7 is, by hypothesis, a simple analytic 
polynomial curve and the Schwarz function S must be analytic on 7 and on its neighborhood (see 
[lUdH]), what exclude the case of the branch point touching the curve 7, concluding the proof. □ 

6 Conformal Deformation 

In this section we conclude the proof of our main result, Theorem [TJ Since the Riemann map h(z; s) 
is conformal from the exterior of the unit disk onto the exterior of j(s) it has a well defined inverse 
from the exterior of 7(5) onto the exterior of the unit disk. We shall denote its inverse by H (z; s): 



for all z in the exterior of 7(5) and s G (0, 1]. 

The Schwarz functions S can be related to the Riemann map h and its inverse H by the 
following: 

Proposition 5 Let 7 be a polynomial curve parametrized by h. Then the Schwarz function is a 
biholomorphic map in a neighborhood 0/7 and is given by 




we would then conclude, contrarily to the hypothesis, that f D f(z)dfi = 0. Therefore, we must 



h(H(z; s);s) = z 



(47) 
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where h(w) = rw + a + a x u> 1 + • • • + a n w n . 
The proof can be found in Refs. [HI ITS] . 



Proposition 6 If h(z; s) satisfies the conditions (i — Hi) of Remark^ then the inverse function 
H(z; s) of h(z; s) reads 

z + \l ' z 2 — AraAs) , . , , . 

H{z- s) = ^ + o(s) (48) 

for < s < e and e small enough. 
Proof. We have by Eqs. (J22D and (J47J) 

n 

h(H(z);s) =rH{z-s) + ^a j (s)H{z-s)~ j = z 

3=0 

which yields 

rH 2 (z; s)-(z + o(s))H(z; s) + o x (s) = . (49) 

Solving the equation for H the result follows. Note that \H\ = 0(1) in the neighborhood of 
7 (s). □ 



Remark 9 Equation (48) has branches, namely the minus square root, which can be directly seen 
from Eq. HOf) . We do not consider the minus square root, because it does not yield S(z) = z on 
7(s). 

Lemma 1 Let dp,(z; s) be the balayage measure supported on 7(5) of the Elbau-Felder ensemble. 
Assume the conditions (i-iii) to hold. Then, given e sufficiently small, for < s < e and f G 
TC(D(s)), we have 

2r 



f(z)dfl(z;s) = —— f{z)y/4r 2 -x 2 dx + o(l) . (50) 



7(s) 2nr 2 J _ 2r 

Proof. Given z in a neighborhood of 7(5) and e > sufficiently small, the Schwarz function of 
7 (s) for < s < e by Proposition [5] reads 

S & s ) = r tj} \ + ai(s)H{z; s) + . . . + a n (s)H n (z; s) 
H(z; s) 

which together with fj48|) . yields 

S(z; s) = Ez + A^z 2 - 4mi(s) + g(z; s) (51) 
where E = (r 2 + \af(s)\)/2rai(s), A = (|a 2 (s)| — r 2 )/2ra 1 (s), and g(z; s) is o(s). By (SSD, we have 

100 = / / m^^dz, 
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for / G TC(D(s)), which together with (l5Tj) . gives 

7 ( s )=/ 7/^( Ez + A v^^M^ + 9(z;s))dz. (52) 

J l{s) 2mt {s) \ J 

The linear term in (1521) will give no contribution by Cauchy theorem. Now we need to estimate 
term depending on g. Since g is analytic in a neighborhood of 7(5), we may have the bound 



27lto{S 



g(z; s)dz 

7 (s) 



< I dz . (53) 



Since arc length §,.dz is finite, max zg7 ( s )|(7(z; s)| = o(s), and by (IToT) and conditions (ii) of 
Remark [3j 

t (s) = r 2 - K(s)| 2 + o(s 2 ) = s + o(s), 

we have that the r.h.s. of ([53]) is o(s)/t (s) = o(l). Hence, A/t (s) = l/2r 2 + O(s). Also, for 

z6 7(s), 

^/4rai(s) - z 2 V4r 2 - ^ 2 
= 2 + ° s • 

Combining the estimates we have Therefore, 

-1 / r . N V4?T 



I ( s ) = TZf /(^ & + o(l) (54) 

«/ 7(5) ' 

By using Proposition ([2]), 7(5) may be deformed until it coincides with the branch cut of the square 
root. Because the branch the integrand along the around the cut will not be the same. In the first 
term running from 2r to — 2r, the square root becomes y/Ar 2 — x 2 , while in the second running 
from — 2r to 2r the square root becomes — -y/4r 2 — x 2 . This leads to 

1 r 2r 

I(s) = / fixWAr 2 - x 2 dx + o(l) 

2?rr 2 J_ 2r 

concluding the proof. □ 



6.1 Proof of Theorem Q} 

We have constructed a one parameter family of curves 7(3) parametrized by h(w; s) with \w\ = 1. 
We have shown, in Proposition [1] that under the hypothesis (1231 the family 7(5) is composed of 
simple closed analytic polynomial curves of degree n. Moreover, h(w; s) acts as the Riemann map 
from the exterior of the unity circle onto the exterior of 7(5) for s G (0, 1]. By Proposition [31 the 
problem of determining the exterior moments tj out of simple closed analytic polynomial curves 
has a unique solution also when t 2 tends to 1/2. 

The balayage techniques enables us to sweep all the eigenvalues to the boundary j(s) of the 
support D(s) and analyze the deformation only focusing on ■y(s). In Theorem 0] we deduce an 
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explicit equation relating the balayage measure with the potential V for / G H. It turns out that 
the balayage measure is proportional to the Schwarz function. Our focus changes to the behavior 
of the Schwarz functions for polynomial curves. 

In Proposition H] we show that the Schwarz function of every simple closed polynomial curve 
must have a branch cut which lies in the interior of the curve and the branch points never touch 
the curve itself. The limit s — > 0, the curve 7 is smashed into the line, see Proposition [TJ The same 
limit is analyzed in Lemma H] where we show that the Schwarz function converges to a function of 
the type \/r 2 — x 2 . 

This implies that as s — > the normal ensemble converges to a Hermitian one (all the eigenvalues 
are real). Using these results and comparing the left hand side of (1501) with the semicircle law (Eq. 
12), we conclude that conformal deformations of Elbau-Felder ensembles to the real line yields a 
Wigner ensemble, as claimed. □ 

7 Examples 

7.1 Potential of degree 2 

Taking n = 1 in <$Ef), the potential V reads 



Consider the one parameter family for Riemann maps h(w; s) = rw+ai(s)w~ 1 in which, without 
loss of generality, r and a^(s) = r(l — s) for s G (0, 1]. Note that this parametrization fulfills the 
conditions (i — Hi) by taking s 1— > s/2. This yields h(w; s) = r(w + (1 — s)w~ 1 ) t = r 2 (2s — s 2 ) 
and 2t 2 = 1 — s, by relations ([IS]) and (ITT)) . The support of the eigenvalues is then given by 



with major and minor semi-axis given, respectively, by (2 — s) and s. For every s G (0, 1] the 
Elbau-Felder conditions is satisfied, namely, t%(s) = and ^(-s)] < 1/2. Condition (|2"5)) : £(s) = 
r — <i2 (s) = s > for s G (0, 1], holds and conditions (i Hi) are trivially satisfied. 

The equilibrium measure associated to this potential is uniform within the support D(s) and 
may be continuous deformed up to the limit lim s ^ h(w; s) = 29fte(u>) for \w\ = 1, in which its 
boundary becomes the line segment [— 2r, 2r] . 

To show that the balayage measure converges to the Wigner measure, we consider the Schwarz 
function, which can be computed using the formulae in Ref. [18] or using H(w; s) together with 
Proposition By equation (1461) and the Cauchy theorem, only the non-holomorphic part of the 



V(z) 



zz — t 2 z 2 — t 2 Z 
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Schwarz function gives a contribution to the integrals. For test functions / G 7i(D(s)) and the 
properties of the balayage measure (|2"Tj) . we have 

/ f(z)dn(z\s) = —yj- r / f(x)y/4r 2 (l-s)-x 2 dx 

JD{s) 2 ^ r (1 - S ) J-2rVT=i 

It is also important to mention that during the whole process of conformal deformation of the 
normal ensemble to the Hermitian, the logarithm potential along the curve is kept unchanged, being 
an invariant in the process of deformation. On can see this by computing the integral explicitly. 



7.2 Potentials of Degree 3: Breaking Down the Hypotheses 

If we break the conditions (i — Hi) of Remark [3] (Aj > 1), then it is impossible to deform the 
normal ensemble into a Hermitian one keeping, at same time, the regularity of the curve j(s) and 
a non trivial support. We shall illustrate this scenario for n = 2 taking t% = = 0. The potential 
then reads 



By (j!6p we have t (s) = r 2 — 2|ci2(s)| 2 and 3ts(s) = a 2 (s)/r 2 . The Riemann map is written as 
h(w;s) = rw + a2(s)w~ 2 . Note that a± = violating condition (ii). Therefore, the deformation 
of 7(s) to a segment of the real line is impossible. The regularity of j(s), which is guaranteed by 
condition (1231) . requires |a 2 | < r/2, while the positive area condition on t requires \a 2 \ < r/y/2. 

There is no parametrization that keeps the regularity of 7. When the area nto of D(s) converges 
to zero (d2 — > r/2), 7(3) develops cusps and is no longer regular. Keeping the regularity of the 
curve, the limit t —>■ is possible only taking r — > 0. Therefore, as the area converges to zero 7 
collapses to a point. 



8 Conclusion 

We study the conformal deformations of Elbau-Felder ensembles to Hermitian ensembles. Special 
attention is paid to Wigner ensembles, for which the density of eigenvalues follows the semicircular 
law. The result presented in Theorem [1] is in a way universal, in the sense that it does not depend 
on the initial ensemble one starts with. Two ingredients are used to prove the main result. We use 
bifurcation theory from a single eigenvalue for extending the moment problem to near slit domains 
and the concept of balayage measure and Schwarz function for assisting the convergence of normal 
ensembles to Wigner ensembles. 

It would be interesting to analyze to deformation of the Elbau-Felder ensembles to the real line 
without the assumption on the convergence rate. In this case the resulting measure, if well defined, 
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might deviate from the semicircular law, and it might also depend on the initial ensemble under 
consideration. 

We would like to acknowledge the financial of the Brazilian agencies FAPESP under grant 
07/04579-2 (TP) and CNPq under 151976/2007-4 (AMV). DHUM is partially supported by CNPq. 



A Solution of (SI) 



Lemma 2 The linear system of equations l[31\) has a unique solution cp = Tv + Bv where B 
(1 — \k\^)~ l JK with k = Kn and T = J + kB, provided \k\ ^ 1. 

Proof. Equation ( 1311 is solvable if and only if 

B T \ ( -K J- 1 \ _ ( I 
f B J y J 1 -K J ~ I / 

holds for some (n + 1) x (n + 1) complex matrices B and T and this is equivalent to 

-BK + TJ- 1 = I 

BJ- 1 -TK = . (55) 

Let us assume that B has its first column given by b and everywhere else (so B has the same form 
of K) and let T = J + kB. Note that KA = a\\K [BA = auB) holds for any matrix A = [a^]. 
Substituting T in ( 1551) . we have 

BJ' 1 -TK = BJ' 1 - (J + kB)K 
= B - JK -\k\ 2 B = 



which implies 



B = l -^JK 

1 - \k\ 2 



I = J h — y IK 



1 - \k\ 

and concludes the proof of the lemma. The uniqueness follows by linearity. □ 
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